

\ 

t 




DEPARTMENT OF MECHANICAL ENGINEERING AND MECHANICS 
SCHOOL OF ENGINEERING , 

OLD DOMINION UNIVERSITY L/J /' 

NORFOLK, VIRGINIA 7 



SURFACE-SLIP EQUATIONS FOR MULTICOMPONENT, 
NONEQUILIBRIUM AIR FLOW 


///'-? 




By 

Roop N. Gupta, 

Carl D. Scott, 
and 

James N. Moss 

Gene Goglia, Principal Investigator 
Final Report 

For the period ending September 30, 1984 



Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 



Under 

Research Grant NAG-1-346 
Ernest V. Zoby, Technical Monitor 
SSD-Aerothermodynamics Branch 


GiLte 


April 1985 


(KASA-CB-181252) SUBFACE-SLJE FCUATIONS FOB 
FC11I CCtf ECNENT, KCtiEgOIlIFBllif AIE FLOS 
Final Report, period ending 3C Sep. 1984 
(Cld Ccicinion Uriv.) 97 p Avail: NfIS HC 
AC5/MF A0 1 CSC! 20D G3/34 


N87-27S45 

Unclas 

C093174 


DEPARTMENT OF MECHANICAL ENGINEERING AND MECHANICS 
SCHOOL OF ENGINEERING 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 


SURFACE-SLIP EQUATIONS FOR MULTICOMPONENT, 
NONEQUILIBRIUM AIR FLOW 


By 

Roop N. Gupta, 

Carl D. Scott, 
and 

James N. Moss 

Gene Goglia, Principal Investigator 


Final Report 

For the period ending September 30, 1984 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


Under 

Research Grant NAG-1 -246 
Ernest V. Zoby, Technical Monitor 
SSD-Aerothermodynamics Branch 




Submitted by the 

Old Dominion University Research Foundation 
P. 0. Box 6369 
Norfolk, Vi rgini a 23508 


April 1985 



TABLE OF CONTENTS 


Page 

NOMENCLATURE iv 

SUMMARY 1 

INTRODUCTION 1 

ANALYSIS 3 

VARIOUS FLUXES AND THE DISTRIBUTION FUNCTIONS 5 

THE BALANCE AND SLIP EQUATIONS 9 

SIMPLIFICATIONS FOR A MULTICOMPONENT MIXTURE 19 

CONCENTRATION SLIP BOUNDARY CONDITION FOR A FULLY CATALYTIC 

AND A NON-CATALYTIC SURFACE . 23 

SIMPLIFICATIONS FOR A BINARY MIXTURE 25 

SIMPLIFICATIONS FOR A SINGLE SPECIES MIXTURE 32 

NO-SLIP BOUNDARY CONDITIONS FOR THE SPECIES CONCENTRATION 34 

DISCUSSION AND CONCLUSIONS 38 

REFERENCES 40 

APPENDIX A: CONSTANTS APPEARING IN THE DISTRIBUTION FUNCTION 

AND DIFFUSION VELOCITY IN TERMS OF TRANSPORT 
PROPERTIES 41 

APPENDIX B: EVALUATED EXPRESSIONS FOR THE FLUXES OF MASS, 

MOMENTUM, AND ENERGY FROM INTEGRALS OVER THE 

DISTRIBUTION FUNCTION 45 

APPENDIX C: SLIP BOUNDARY CONDITIONS GIVEN BY DAVIS FOR A 

BINARY MIXTURE 54 

APPENDIX D: SLIP CONDITIONS IN THE BODY-ORIENTED COORDINATE 

SYSTEM 57 

APPENDIX E: SLIP CONDITIONS IN THE SPHERICAL COORDINATES 79 

APPENDIX F: INTEGRALS REQUIRED FOR OBTAINING VARIOUS FLUXES 85 

ii 


TABLE OF CONTENTS - Concluded 


LIST OF FIGURES 

Figure Page 

1 The Knudsen layer showing general fluxes and coordinate 

axes. The temperature as a function of normal distance 
is schematically overlayed 36 

D1 Coordinate configuration for a Body (Symmetric about its 

Axis) 57 

El Coordinate configuration 79 


iii 


NOMENCLATURE 


A 1 R 1 r . 1 ( J ) 

A k » B ki * L k 

a iO’ a il» b iO’ C S^ 


Coefficients of the velocity distribution 
function for ith species 

Distribution function coefficients defined in 
ref. 5 for ith species 

Speed of sound, 7 R T « 

Mass fraction of species i 


C pi Specific heat at constant pressure of species i 

Cp Frozen spej^fic heat at constant pressure of the 

mixture, ICC 
i i i 

D.,- .• Multicomponent diffusion coefficient for species 

pair "i" and "j" 


^ij 


-i 


E' 

e 


F 

f 1 

f(0)1 


h 

h 

K 

k 


Binary diffusion coefficient 

Thermal diffusion coefficient 

Diffusion vector of .ith species 

Total energy flux, E T + eM 

Energy flux from translational energy 

Internal energy that readily equilibrates with 
translational energy 

General flux of property such as mass, momentum, 
or energy 

Distribution function of ith species 
Equilibrium (or Maxwellian) distribution function 
of ith species 
NS 

Enthalpy, z h C 
i=l i i 

Enthalpy of species i 
Thermal conductivity of mixture 


Boltzmann constant 



'W1 

L ij 

£ 1j 

Le 

M 

m i 

m 

NS 

"i 


n 


Pr 

P 

P 

Pi 

R 

★ 

R 

Re 

R 

r 

r N 

* 

s 


Wall catalytic recombination rate constant for 
species i 

Multicomponent Lewis number, Cp p D^j/K 
Binary Lewis number, C p p ^ j/K 
Binary Lewis number, C p p D m /K 
Mass flux 

Mach number, U^/a*, 

Mass of ith species 
Mass of a mixture molecule 
Number of chemical species 
Number density of species i 
Total number density 

Coordinate measured normal to the body in the 
body-filled coordinate system (Appendices C 
and D) 

Prandtl number, CpU/K 
Momentum flux 
Pressure 

Partial pressure of species i 
Wall heat-transfer rate 
Gas constant 

Radius measured from axis of symmetry to a point 
on the body surface 

Reynolds number, U. r N /u pef 

Universal gas constant 

Radial distance from body center (Appendix E) 

Body nose radius of curvature (Appendix D) or 
radius of sphere (Appendix E) 

Distance measured along the body surface 
(Appendix D) 


v 



T 

U„ 

u 


vj 

K 

v ok 

w ok 

w k 

w i 

R 

x 

y 

z 

a 

8 

5 k£ 

e 

Y i 

7 

K 

\ 

V 

$ 

♦ 1 (^) 


Temperature 

Free stream velocity 

Velocity component along the body surface 
(Appendices C and D) 

Thermal (or peculiar, or molecular) velocity of 
the ith species 

Diffusion velocity 

Bulk or mass averaged velocity 


Nondimensional mass averaged velocity, 
Nondimensional thermal veloci 


& 
V 2KT 


ity ’ vir 

Molecular (or atomic) weight of species i 
Molecular weight of mixture 
Coordinate parallel to body (Fig. 1) 
Coordinate perpendicular to body (Fig. i). 
Coordinate parallel to body (Fig. 1) 

Energy accommodation coefficient 
Body angle defined in Fig. D1 
Kronecker delta 

i 1 2 

Reynolds number parameter, [u re f/p 00 U 00 r N ] 1 

Recombination coefficient (or, fraction of 
incident atoms that recombine at surface) 

Ratio of specific heats 

Surface curvature 

Mean free path 

Viscosity of mixture 

Perturbation part of distribution function 


'ok 


General property of 
momentum and energy 


ith species such as mass. 



<{> 

p 

9 » 9 i 


Circumferential angle (Appendices 0 and E) 
Density 

Accommodation coefficient (or fraction of 
incident particles that stick) 


Affixes (Used both as subscripts and superscripts) 


i» J 

s 

w 

Subscripts 


Species indices 

Edge of Knudsen layer 

Wall 


A 

e 

M 

i, j 


Atom 

Equilibrium value 
Molecule 
Specie indices 


q 

ref 

x, y, z 

y 

I ! , y, z 


Co 0?’d 1 Hits 1 (J 1 C wS 

Specie index 
Reference quantity 
Component directions 
Normal component 
Tangential components 
Freestream value 


Superscripts 


+ 

+ 


Dimensionless quantity; also used locally 

Incident flux 

Specularly reflected flux 



SURFACE-SLIP EQUATIONS FOR MULTICOMPONENT, 

NONEQUILIBRIUM AIR FLOW* 

Roop N. Gupta** 

Carl D. Scottt 
and 

James N. Mosstt 
SUMMARY 

Equations are presented for the surface slip (or jump) values of spe- 
cies concentration, pressure, velocity, and temperature in the low-Reynolds 
number, high-altitude flight regime of a space vehicle. These are obtained 
from closed-form solutions of the mass, momentum, and energy flux equations 
using the Chapman-Enskog velocity distribution function. This function 
represents a solution of the Boltzmann equation in the Navier-Stokes ap- 
proximation. The analysis, obtained for nonequilibrium multicomponent air 
flow, includes the finite-rate surface catalytic recombination and changes 
in the internal energy during reflection from the surface. Expressions for 
the various slip quantities have been obtained in a form which can readily 
be employed in flow-field computations. A consistent set of equations is 
provided for multicomponent, binary, and single species mixtures. Expres- 
sion is also provided for the finite-rate species-concentration boundary 
condition for a multicomponent mixture in absence of slip. 

INTRODUCTION 

For an accurate prediction of the aerothermal environment of a space 
vehicle entering the earth's atmosphere in the high-altitude low-Reynolds- 
number flight regime (Ref. 1), the multicomponent, nonequilibrium gas 
chemistry as well as the wall slip and catalysis effects, must be evaluated 

*The authors are grateful to Dr. Fred G. Blottner (of Sandia National 
Laboratories, Albuquerque, NM) for many useful discussions pertaining to 
the mass flux expressions contained in Appendix B. 

**01 d Dominion University, Norfolk, VA 
t Johnson Space Center, Houston, Texas 
ttLangley Research Center, Hampton, VA. 



Such flow fields are of particular interest for aeroassist and space shuttle 
vehicles. This study was undertaken to obtain the boundary relations that 
incorporate the effects of slip, multicomponent diffusion, wall catalycity, 
and changes in internal energy of the molecules (during reflection from the 
surface) for application to flow-field calculations under the general as- 
sumption of local thermodynamic equilibrium. 

Scott (Ref. 2) first presented the wall boundary conditions for a 
multicomponent mixture with diffusion and wall-catalysed atom recombination. 
In obtaining these boundary conditions, he used a first order velocity dis- 
tribution function at the edge of the Knudsen layer next to the wall, where 
the continuum model of the gas breaks down. These boundary conditions, 
obtained from the kinetic theory considerations, provide solutions at the 
top of the Knudsen layer that would match the solution of the Navier-Stokes 
equations in the bulk outer flow. Hendricks (Ref. 3), using Scott's formu- 
lation, obtained simplified expressions for engineering applications with 
some corrections to Scott's expressions. Hendricks' analysis, however, 
contained some gross errors in obtaining the engineering expressions. This 
paper reanalyses the wall boundary equations by using the approach of refer- 
ence 2 and provides appropriate relations for the various quantities with 
surface slip in a form which can readily be employed for flow-field computa- 
tions. An effort has also been made to reconcile the differences between 
slip expressions employed by the different researchers. The present analy- 
sis provides a consistent formulation for the slip equations for a multi- 
component, binary, and single species mixture. Expression is also provided 
for the finite-rate species-concentration boundary condition for a multi- 
component mixture (in addition to that for a binary mixture) in the absence 
of slip. This may be of interest for the shuttle flow-field calculations 
based on multicomponent diffusion (Ref. 4). 
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The main difference (among other details) between the results contained 
in reference 5 and the present work is that the various internal degrees of 
freedom for a molecule were considered frozen during reflection from the 
surface in reference 5, whereas, they are allowed to change in the present 
analysis.*** Further, the results are provided in both spherical and body- 
oriented coordinate system. 

ANALYSIS 

The slip conditions are taken to exist across the Knudsen layer, which 
is on the order of one mean free path in thickness as sketched in Figure 1. 
The analysis outlined here follows the approach of References 2 and 5. It 
is based on Shidlovskiy' s (Ref. 6) assumption that the distribution function 
near the wall can be described to first-order accuracy by the so-called 
Navier-Stokes approximation. However, a deviation is made from the proce- 
dures of Shidlovskiy in that a Chapman-Enskog type distribution function for 
a multicomponent mixture obtained by the variational method of Hi rschfelder, 
Curtiss, and Bird (Ref. 7) is used. The Chapman-Enskog distribution 
function allows accounting for diffusion. The analysis contains the follow- 
ing assumptions: 

(i) The energy and momentum accommodation coefficients (i.e. a 
and e, respectively) have the same value. 

(ii) The fluxes of mass, momentum, and energy across the Knudsen 
layer are assumed constant. This is consistent with the 
assumption of negligible variation of the velocity distri- 
bution function through the Knudsen layer. 

***Since the assumption of local thermodynamic equilibrium is employed, only 
those internal energies are considered which equilibrate readily with the 
translational energy. 
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(iii) The internal energy associated with the rotational and 

vibrational modes readily equil ibrates with translational 
energy. 

The interaction model at the gas-solid interface, with the various 
fluxes sketched in Figure 1, can be mathematically stated for dissociated 
air as provided through the following equations. 

For a recombining atom : 

F A = F A + F A + ■( 0 A‘ Y A ) F A ; A = 0 ’ N (D 


For a molecule gaining from the corresponding atom recombination : 


F M = F M + + 9 m F m + Y A F fl ; M = 0 2 , N 2 


M ' M M M 'A A 


( 2 ) 


* For all other atoms and molecules (surface is assumed to be noncata- 
lytic with respect to them): 


F, ' A * < lJS 1> F * + e iF 1 " 


(3) 


where F^ denotes a convective property such as mass, momentum, or energy. 

Summing over all the species, we obtain from Eq. (1), (2), and (3), the 
following expression for the net flux of moment in or energy: 


NS 


l 

i=l 


F i 


NS , NS 

I F + . + I (1-6 ) 
i=l 1 i=l 1 


F+ i 


NS 


+ I 

i=l 


9 i F i 


(4) 



That is, the net flux at the outer edge of the Knudsen layer equals the 
incident flux, plus the specularly reflected flux (incident minus the frac- 
tion that sticks) at the wall, plus the diffusely reflected flux (those that 
accommodate to the wall) from the wall. 

Each species is treated separately in the mass balance equations. 
Therefore, Eqs. (1), (2), or (3) are employed depending on the species being 
considered. In Eq. (1), the diffusely reflected flux consists of those 
atoms that are accommodated to the wall minus those that recombine. For the 
molecules in Eq. (2), the diffusely reflected term is present along with the 
source term resulting from the appropriate atoms recombining on the surface. 

VARIOUS FLUXES AND THE DISTRIBUTION FUNCTIONS 
The interaction model of Eqs. (1) through (4) is employed to obtain the 
slip boundary equations at the gas/solid interface. Through Eq.(l), the net 
fluxes of species, momentum, and energy at the outer edge of the Knudsen 
layer are equated to the difference between the incident and reflected 
fluxes at the wall. These fluxes are assumed to be constant across the 
Knudsen layer and are obtained from moments of the distribution function. 

For a convected property ^(v) such as mass, momentum, or energy for the 
ith species, the net flux of that property normal to the wall at the outer 
edge of the Knudsen layer, for example, is 

Fi «L" ' C C VyOVv) f s (*> d3vi (5) 

where is the normal component of the molecular velocity and f^ is the 
velocity distribution function at the edge of the Knudsen layer. 

Similar integrals are obtained for the incident and reflected fluxes by 
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integrations over appropriate half-spaces in molecular velocity: 


Incident flux: 


Fj ■£ L C *',*'<*) f’(5) dV (6) 

Specularly reflected flux: 

’ C Q C Vy +*(5) fs(V x . -* y . VyJd^V 1 (7) 

Diffusely reflected flux: 

F"=£ r 0 £ v’^di) f’(5) d^v 1 (8) 

where f^ is the Maxwellian velocity distribution function. 

The velocity distribution functions used in integrals contained in 
relations (5) through (8) are those for a nonuniform multicomponent mixture 
perturbed out of equilibrium: 

fi(^) = f W'fi) [1 + ^ (V)] (9) 

where f^^(V) is the Maxwellian distribution function for the ith species 
given as 


f (0)i (1 1) = 


n i m i 


{2irkT) 


3/2 

3/2 


-(m./2kT)V. 2 

e 


( 10 ) 
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and 


♦Yv) = - 


3 V 


- B 


Ok 


3 X, 


kz 


3x. 


NS 
+ n I 
j=l 


•i( j) 


i 


(id 


i 


Here k and Z are the dummy indices for three coordinate directions and 
the summations with index j represents summation over all the species. 
The summation convention for repeated indices is used. The coefficients 
A|<, anc * ^k^ are functions of the dimensionless velocity: 



( 12 ) 


a nr! :ro Hafinpjj 

wkt ■ u vk( w uw I i ilQv u J 



(13) 



(14) 



_(j) U 1 

c i0 k 


(15) 


where a^, a.p b^g, and c^g ; are constants determined from the varia- 
tional problem in the first approximation for a mixture as given in ref. 7 
and S|^ is the Kronecker delta such that 


<5. . = 1, if k = * 
k z 

6 k* = °> if k * l - 
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These constants are functions of the collision integrals and are related to 
the transport properties. The simplified form of these constants is pro- 
vided in Appendix A of this paper. More details can be found in Appendix A 
of Ref. 2, or in Ref. 7. The diffusion vector d^ in Eq. (11) is related 
to the diffusion velocity of the jth species and is defined after neglecting 
the external forces as (Ref. 7): 


d 


j 

k 


. n . n . n . m . » 

— ( J.) «■ ( J- - -M — (»np) 
3 x k n ii p 3 x k 


(16) 


where n . and m. are, respectively, the number density and mass of the 

J J 

jth species, n is the total nunber density n = £ n-, p is the total 

j J 

mass density p = I n.m. , and p is the total pressure p = £ p. . A 

*4 J J [ 

-• . j 

simplified form for d k is provided in Appendix A. 

The total mass averaged velocity v Qk (i.e. the kth component) 
appearing in Eq. (11) is defined as 


Ok 


ns 

l 

j=l 


n .m . 
J J 
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where v k is the mean of total velocity v k = Vg k + of the jth 
species averaged over the distribution function and V k is the thermal (or 
peculiar) velocity, also introduced in Eq. (5). The thermal velocity V k 
of the jth species averaged over the distribution function is known as the 
diffusion velocity V^: 
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C C C vj (ty f j (^) d 3 v j 


(17) 


vj = — 

"j 


A simplified expression for the diffusion velocity in terms of transport 
properties is provided in Appendix A. 

THE BALANCE AND SLIP EQUATIONS 

Based on Eqs. (1) through (4), the balance equations for the ith 
species for fluxes normal to the surface of the species mass the normal 
component of momentum P^, the tangential component of momentum P-j||> an d 
the energy E^ are obtained as follows: 

(a) Species mass flux 

For a recombining atom 

M Ay s 9 A M A + (0 A _ 7 a) ^ ? A = 0, N (18) 

For a molecule gaining from the corresponding atom recombination 

X ' 9 M < * S M ^ + *A*A i N 2 (19) 

For all other atoms and molecules 



* 9 i H i + e i< 


( 20 ) 


(b) Normal momentum flux 



( 21 ) 
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(c) Tangential momentum flux 


NS 

I P, 

i=l 


i | 


NS 


l 9 . P. 


i=l 


l 1 


( 22 ) 


(d) Energy flux' 


NS 

I E< 


i=l 




NS 

I 

1-1 


e t + 
iy 


I 


e. M. 

iy iy 


NS Tl ns 

I Mi** I 

1=1 1 1-1 


diatomic 

molecules 

Tw 


8 lV" ♦ l e.eft 


1 

diatomic 

molecules 

„ w„w 
6, e, M, 


diatomic 

molecules 


(23) 


where e^ is the internal energy of ith species that readily equilibrates 
with the translational energy e| under the assumption of local thermo- 
dynamic equilibrium. For example, 

. rotation kT 


In writing the energy flux balance of Eq. (23) it has been assumed that 
there is no change in internal energy during specular reflection. 

In obtaining Eqs. (18) through (23) we have used the following 
rel ations 


tThe energy balance is based on the assumption that the various energies 
considered readily equilibrate with the translational energy. 
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3 




and 



Because it is assumed that the atoms are consumed at the wall by catalytic 
recombination in Eq. (18), the net mass flux M^O. Similarly, the net 
mass flux M^y * 0 in Eq. (19). However, M. = 0 in Eq. (20) for the 
atoms and molecules for whom the surface is assumed noncatalytic . 

By substituting Eqs. (5), (6), (7) and (8) with the definition of 
♦ 1 ( V) as mass for Eqs. (18) through (20), as normal component of momentum 
for Eq. (21), as tangential component of momentum for Eq. (22) and as energy 
for Eq. (23), respectively, and carrying out the integrations*, one obtains 
equations relating the slip properties to wall properties and gradients at 
the edge of Knudsen layer. All accommodation coefficients 0^ are assumed 
to be equal to 0. 

Number density (or concentration) slip (obtained from mass flux balance): 



Pressure slip (obtained from the flux balance of normal component of momen- 
tum) 


tThe mass, momentum, and energy fluxes in terms of evaluated integrals over 
the distribution function are given in Appendix B. 
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„ 3 3 V 3 V 

p s - (-- [» (_ 2 i_ + . 2 _oy)j 

3 3 x 3 z 3y S 


+ P- [k- f n,(a 10 -la )] 
W 3y i = l 1 1U 2 11 


+ i p w 


/I - *- I [n, I 

2 VtT i 1 -- 


NS 


9 . 2-0 *■ r_ r~’ (j) .Ji i 

j-1 10 ^ 


( 25 ) 


. NS 

with viscosity n = i kT l n- b, rt 

2 i=l 1 10 


Velocity slip (obtained from the flux balance of tangential component of 
momentum) 


S _{-/-( 2-0'i r V 0x V 0y-, 
v 0 = { w (_j u ( + — q 

20 s 3y 3 x s 


NS 


*- i p* [ ®£nl 

2 1-1 1 3 x 


U i0 - J a il) - " c So’ d xls> 


NS m. 

/ l (Pf -1~) 
i=l 1 J 2k T 


(26) 


V0 Z = {^ . (— ) M.( 


_ 3 v 3 v 

2-0i .. r Oz . Oy^ 


20 


3y 3 z 


, NS _ , NS ... . 

+ ? ! 7T (a 10 '7 ai 1 ) '-/^c'o d z ] s> 
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Temperature slip: 



Equations (24) through (28) differ from the corresponding expressions 
provided in Ref. 2 due to small errors and also due to the differences in 
the interaction model employed at the gas-solid interface. 

The constants a i0 , a.^, b^g, and c^ (also known as the Sonine ex- 
pansion coefficients) appearing in Eqs. (24) through (28) may now be ex- 
pressed in terms of the transport properties as given in Appendix A. Using 
these relations along with the various flux expressions of Appendix B and 
also expressing d^ in terms of the gradient of mass fractions 3C./3y as 
given in Appendix A (by neglecting diffusion due to pressure gradients and 
external forces), the following equations are obtained after some algebraic 
simpl ifications:^ 


§The approximations made in the expressions for a.-i and bin are given in 
Appendix A. 11 10 



Number density (or concentration) slip 


Same as Eq. (24). 


Pressure slip 


P s .{.lu + !^. 2 !>)] s + i (i±) [I »i f „. Kl E 

3 3x 3 z 3y s 5 YT = 3y i*l 1 1 yj 2kT 
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Velocity Slip 
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NS ■— _ NS i — NS 

h "1 h K i ] s - S s ^ ®u [ 

j*i 


!i 

3 Z 


NS S 3C „ NS _ 

- Cj I ^ _i)] } / X nf 

J q=l m„ 3 z i=l 1 \ 1 


( 31 ) 


Temperature slip 
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( 32 ) 
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where the mass fraction and the mass of a mixture molecule m are 
defined, respectively, as 


and 

NS _i 

m = ( I C ./m.) 
j=l 3 3 

Equations (24), (29), (30), (31), and (32) differ from Eqs. (25) through 
(28) of Reference 3. For dissociated air, the gas-surface interaction model 
employed in Reference 3 appears to be inconsistent. Further, the diffusion 
vector d^ appears to have been incorrectly evaluated in Reference 3. If 
one carries out the simplifications in Eqs. (23), (24), and (28) of Refer- 
ence 3 (which contain the binary assumption) through the evaluation of 

g 

— (n./n) in terms of mass fractions as given by Eqs. (35), (36), and (39) 

3y J 

of the same reference, erroneous results are obviously obtained because 
NS 

I (3C-/3r) (which is zero by definition) is contained as a factor in 
j=l J * 

several of these equations. 

For the first-order recombination at the surface, the following 
relation between the atom mass flux and the wall number density n^ 

may be employed 


M Ay = " k wA ^ n A V ( 33) 

where the minus sign indicates that the flux is in the direction opposite to 
the outward normal and the expression for the rate constant k wA with 
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diffusion and slip is (Ref. 2). 


tt 



( 34 ) 


Here y ^ is the recombination coefficient. 

For a fully catalytic wall (y^= 1), the maximum value of the rate 
constant k . is limited by the surface temperature. The reaction rate con- 
stant for a fully catalytic wall with the gas phase in chemical equilibrium 
is often assumed to be infinity for the sake of simplicity. 

Employing Eq. (33), the net mass flux M.. appearing in Eqs. (24) and 
(32) may be defined as 



“ k wA n A m A ; 


A = 0, N 




„w m . 

A V 


r A = 0 for M = O 2 
A = N for M = N 2 


For all other atoms and molecules. 


(35a) 


(35b) 



(35c) 


Equation (24), with M. defined by Eqs. (35a) through (35c), gives 


TtAs shown in ref. 2, by neglecting slip but keeping diffusion, a slightly 
different form of Eq. (34) is obtained: 


2 T a 

k wA - (— ) 

2-V, 


kT. 


w 


2irm n 
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the number density ratio n^/n^. However, to obtain nV from this ratio, 
an expression for n* is required. 

The net mass flux of 0 and N atoms to the surface, M^, is also 
equal to the rate of consumption of these atoms at the wall from surface 
recombination: 


M Ay = y A% S A = °* N 


(36a) 


and the corresponding net mass flux of O 2 and N 2 molecules will be 


M = - V M + • |A = 0 for M = 0 2 

My y A n Ay »’ ‘A - N for M » N 2 


(36b) 


For all other species, the net mass flux to the surface may be assumed to be 
zero. Thus, 

M 1y = 0 (36c) 

Substituting values of the net and incident fluxes M.^, M.. from Appendix 

8 in Eq. (36), the following expression is obtained (after we neglect 
thermal diffusion): 


J : 


NS e 9 C . NS = 9 C 

I [— - c- l (5 !)],*- 

3=1 J 9y J q=l m„ 9y m s 


(37) 
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which may be used to obtain 



Here ip' 


is the source term defined as: 


,s r A 
♦a -(Z=r}J 


i m.n. s / 2kT p. 

I JJL / !_ ( _ii + l) 

VIT 2S s V m A Pft 


(38a) 


For Q 2 and N 2 molecules 



r A = 0 for M = O 2 
1 A = N for M = N 2 


(38b) 


For all other species 

+ • = 0 


(38c) 


SIMPLIFICATIONS FOR A MULTICOMPONENT MIXTURE 
Equations (24) (29), (30), (31), and (32) for multicomponent gas flows 
can be simplified if one makes the following assumptions: 


(i) All the diffusion coefficients, D.., for a multicomponent gas have 

* J 

§If no assumption is made about Da,-, Eq. (37) would give an expression for 
(3C i /3y) $ for all the species: J 


3C. m. 

Cr^s-C^ { 

3y * m,. 


NS 
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4=1 
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'ij 
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L- 1 - c j 
3y J 


NS £ 
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q=l m c 


3C NS 
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ij c j)s 


The source term ip.j in this expression, however, may be simplified to yield 
an expression for C| (or n|) for the recombining atoms only : 


r s m A n A 
U A " 


2-y r m A 

- 2 VT (_ 1 ) (J.) 

y A m_ 
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ft NS 0 3 C • NS = 3C 
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is the same as the 


the same value so that D. . = D 12 ,* D 12 

J 

binary diffusion coefficient D 12 . 

(ii) The normal momentum flux to the pressure ratio, P-jy/Pj is the 
same for all species and equal to that of the mixture. This also 
implies that the normal shear stress t . for species i is the 
same as that for the mixture (t ). 

v yy' 

(iii) The rotational and vibrational states are fully excited so that 
the internal energy e^ for the air molecule may be taken as 
equal to 2kT/m.... 


*This is a somewhat stronger assumption. Because D-j-j are concentration de- 
pendent, whereas D-j,- are virtually independent of composition. The multi- 
component diffusion J coeff icient • is related to the binary diffusion 

coefficient Z)-j^ through the following relations (see Ref. 7): 


D ij ■ K ij - (W *ii 


where quantities K<j are coefficients in a matrix which is the inverse of 
the matrix with the following coefficients: 



A 

^ij 

0 


NS 

M 0 Z. 

J q=l 
q^i 



(i*j) 

(i=j) 


One can see, therefore, that by employing D 12 (which is same as ^ 12 ) 
for all the species in a multicomponent gas mixture, considerable saving is 
obtained in computational effort and time without losing the general flavor 
of multicomponent diffusion. This is particularly true if the dissociated 
air consists predominantly of nitrogen molecules and oxygen atoms. The 
shuttle entry conditions fall into this category. Simple and multicompo- 
nent diffusion gave same results in "An Experimental and Analytical Study 
of Slip and Catalytic Boundary Conditions Applied to Spheres in Low 
Reynolds Number Arc Jet Flows," by Carl D. Scott, proceedings of the 9th 
International Symposium on Rarefied Gas Dynamics, Gottingen, July 15-20, 
1974, pp. D.14-1 to 11. 
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These assumptions retain the major effects of multicomponent fluxes on 
various slip quantities and provide considerable saving in computational 
effort required for the analysis of a flow problem (ref. 2 and footnote on 
page 20). With these simplifications, Eqs. (24), (29), (30), (31) and (32) 
yield: 


Concentration slip 
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Pressure slip 
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Velocity slip 
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Temperature slip 
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molecules (43a) 
or, if the internal energy is frozen during reflection from the surface (see 
Eq. (28) of Ref. 5, for example). 
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Equation (37) may also be simplified to yield an explicit expression for n^ 


where 


-I 


ac. , ii. 

(—1) + — (— ) s 
p. ay m. d 

_i [ 1 + ] i£ — ] 

m, NS = 3 C 

l (5 !) s 

q=l m„ 8 y 
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(44) 


*A = - 


1 m A n A 


(2-yJ ^ 


2n. 


/2kT P. 

— (— + 1) ; A = 0, N 

m A Pa 


(45a) 


s _ ^s . { A = 0 for M = 02 

M A * 1 A = N for M = N 2 


(45b) 


= 0 for all other species (45c) 

It is suggested here that the concentration for the major specie (for 
example, nitrogen) be obtained by requiring the sun of concentrations of all 
the species to equal unity. It may be mentioned here that the mass of the 
ith species, m., is related to the molecular (or atomic) weight, VL , 
through the relation 


— =- (46) 

W. R 

where k is the Boltzmann constant and R is the universal gas con- 
stant. 


CONCENTRATION SLIP BOUNDARY CONDITION FOR A FULLY CATALYTIC 
AND A NONCATALYTIC SURFACE 

Equation (44) gives slip values of the concentration n* for a finite 
catalytic surface. For a fully catalytic (y a = 1) surface one generally 
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assumes complete recombination of atoms at the surface. There is a slight 
discrepancy in this assumption because the maximum recombination rate is 
limited by the surface temperature as discussed earlier [See the discussion 
following Eq. (34)]. Thus, for a fully catalytic surface, equation (44) or 
(41) should be employed with y^ = 1. 

For a noncatalytic surface (y^ = 0), Eq. (45) gives = 0 for all 
the species. For this case, Eq. (44) becomes 


s ^ s r i 
* — [1 + 


ac. 

3y 


m i 


NS j 3 C 

I (S- -A 

q=l m q ay 


(47a) 


which may also be written as 


3C. 

(— ). 
ay 


c NS s 3 c 

S\ r rm Ch 


- - (1-c-) l { 


q=l m q 3y 


(47b) 


Summing the above equation over all the species gives 


NS 5 3 C 

I — ) s = 0 

q=l m n 3y 


(48) 


for a noncatalytic surface. Therefore, Eq. (47a) is not an appropriate 
boundary condition. However, employing Eq. (48) with Eq. (47b) yields: 


(— L = 0 (49) 

ay * 
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which may be used as the boundary condition for a noncatalytic surface with 
a multicomponent gas mixture. 

SIMPLIFICATIONS FOR A BINARY MIXTURE 


At lower altitudes with the flow in slight nonequilibrium (i.e. when 
the Reynolds number parameter e 2 << 1), the derivatives of various quantities 
with respect to x and z and some other higher order terms like 8 V 8y 
may be dropped through an order of magnitude analysis (see Ref. 9, for exam- 
ple). Further, the dissociated air may be considered as a binary mixture 
(i.e. consisting of atoms and molecules only, see Ref. 8) at these altitu- 
des. If an assumption is also made that the internal energy of the mole- 
cules remains frozen during reflection from the surface, Eqs. (39) through 
(45) can be simplified further to the forms given here. 
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or, with frozen internal energy during surface reflection. 
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In Eqs. (54a) and (54b), n A is obtained from n A by using Eq. (50). 


Expression for n] 
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which may also be written as 
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(55b) 


n A m A 


= c- f 



s 

^AM 
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For a fully catalytic surface (y^ = 1), Eq. (55b) gives appropriate 
value for the concentration slip.^" The relevant boundary condition for a 
noncatalytic surface (y a = 0) will be 

3C. 

(— L = 0 (56a) 

3y s 


as can readily be seen from Eq. (55b). 

For a noncatalytic surface (y^ = 0), pressure and temperature slip Eqs. 
(51) and (54b) are further simplified, with the help of Eq. (56a), to: 


Pressure slip 
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Temperature slip (with frozen internal energy during reflection from the 
suface) 
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(58a) 


Sometimes in the literature is prescribed as zero for a fully cata- 
lytic surface. Strictly speaking, this will be true only when the 
Reynolds number parameter e 2 (Ref. 8) is approximately zero (close to the 
chemical equilibrium condition at low altitudes). This can easily be seen 
by nondimensionalising equation (55b) in a way similar to equation (28d) of 
reference 8. Thus, the recombination rate coefficient y/\ and density (as 
measured through e 2 ) , both control the recombination rate and not y/\ alone. 
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Equations (52a), (53a), (57a), and (58a) ( with the noncatalytic surface 
assumption in Eqs. (57a) and (58a)) can be rewritten as: 


Velocity slip 


Pressure slip 


Temperature slip 
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a l = 


5rr 

16 


fl (i±) = 1.2304 (H-) 
>26 9 


“ f- ft ■ i. 

16 ^2 9 


1750 (£i) 
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IS! [T (?±) . 2.3071 fill 
128 > 2 9 0 


The concentration slip condition consistent with Eqs. (52b), (53b), 
(57b), and (58b) is 


ac 

L = 0 (56b) 

3y 


The range of values of Ai , Bi , and Cl is as given below for a mixture of 
oxygen atoms and molecules and y = 1.4 (ref. 8): 


1.0039 < Ai < 1.0186 

0.9507 < B 1 < (0.9507 )x(Jl — +■ JT + 

4 m A 2 m M m A m M 

0.9056 < Ci < 0.9507 

The minimum values for Ai and Ci occur at n^ = 0.5 P/m^, whereas the 
maximun value for Bi occurs at this value of n^. 

The expression (52b), (57b), and (58b) reduce to those obtained in Ref. 
8** if one assumes the values for \ , B x , and C L to be unity and a non 

*^See Appendix C for the dimensional form of the slip boundary conditions 
given in Ref. 8. 
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catalytic surface boundary condition as given by Eq. (56b). The concentra- 
tion slip boundary condition provided in Ref. 8 is: 



which is the same as Eq. (55b) obtained here except for the additive term 
s ttt 

C^. It is obvious that the temperature slip boundary condition of Eq. 
(2.8b) obtained in Ref. 8 is valid, strictly speaking, for a noncatalytic 
surface only and is not consistent with the concentration slip boundary 
condition of Eq. (2.8c) (reproduced here as Eq. (59)) obtained for a finite 
catalytic wall. As a matter of fact, the velocity and temperature slip 
boundary conditions of Ref. 8 are similar to those of Ref. 9,^ where these 
slip conditions are provided for a perfect gas (or single species mixture). 

The inconsistencies in the boundary conditions used in Refs. 10 and 11 
are similar to those of Ref. 8, namely, the pressure and temperature slip 


^^It appears that Eq. (2.8d) of Ref. 8 for the concentration slip has been 
formulated for small deviations from the chemical equilibrium condition 
(i.e. the Reynolds number parameter e 2 << 1). Thus, when flow goes to 
chemical equilibrium (with e 2 « 0), one obtains from Eq. (2.8d) the 
equilibrium value for the concentration i.e. which would be 

zero for the oxygen atoms for surface temperatures of 2000°K or less. 

^There appears to be some error with the form of Eqs. (2.7c) and (2.7d) 
given in Ref. 9 if one employs the definition of dimensionless heat-trans- 
fer rate ' q ' in these equations from Eq. (2.4b) of the same reference. 

(The Reynolds number parameter, e 2 , is missing in Eqs. (2.2b) and (2.4b) 
for dimensionless shear stress and heat-transfer, respectively, of the same 
reference. These have been corrected in Ref. 8.) The resulting pressure 
and temperature slip expressions will contain the dimensionless viscosity 
coefficient 'u 1 without the Prandtl number, see Ref. 10 and 11, for exam- 
ple. This is in contradiction to Eq. (2.8c) of Ref. 8, which contains 
dimensionless * K * in place of dimensionless 'y 1 . For a perfect gas, 
employing 'y ' in place of 1 K ' in the pressure and temperature slip equa- 
tions will result in no error. However, for real gas flows, the pressure 
and temperature slip values will be in error by the factor of y/K. The 
present analysis indicates that Eq. (2.8c) of Ref. 8 is of the "correct" 
form for a single species mixture. 
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values are in error by a factor of dimensionless viscosity to the thermal 
conductivity ratio 'u/K'. References 10 and 11 also employ the single 
species formulation for obtaining the velocity, temperature and pressure 
slip values and the concentration slip expression in these references is 
that obtained for a binary mixture (Ref. 8). The coefficients appearing in 

the slip equations of Ref. 11 may be obtained by multiplying ai , bi , and 
q (coefficients of Ref. 8) given here by 16/5-ir. The factor 16/5 is some- 
times repl aced by it . 

SIMPLIFICATIONS FOR A SINGLE SPECIES MIXTURE 

Eqs. (39) through (42) and Eq. (43b) may be simplified for a single 
species mixture to the expressions obtained in Ref. 6. For small jump (or 
slip) conditions these simplified equations may be written as (with Y A = 0 
and M. - 0 for a single species mixture): 

Density slip 
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where we have used the following relations between the coefficients of vis- 
cosity and thermal conductivity and the mean free path for perfectly elastic 
spherical molecules (ref. 12). 


y = -L_ jzn RT pX 

16 V 


(65a) 
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15 k 
— _ y 
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(65b) 


The expression for y given earlier, following Eq. (25), for a single spe- 
cies mixture is modified to: 

y = 1 nkTb 0 
2 


If, now, an assumption is made that the slip values are small , equations for 
the density, pressure, and temperature slip are simplified to: 


* Si nee the temperature slip expression contained in Ref. 6 is for a gas 
consisting of perfectly elastic spherical molecules, the gas does not 
possess any internal degree of freedom. Therefore, we have obtained the 
temperature slip equation from Eq. (43b). 
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Density si ip 
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In obtaining Eq. (64c) from (64b), we have used P y /p s « 1 in the second- 
order terms (i.e. 2nd and 3rd terms on right side of Eq. (64b)). Equations 

(61b), (62), (63), and (64C) are the slip equations given in refs. 6 and 

12 . 

NO-SLIP BOUNDARY CONDITION FOR THE SPECIES CONCENTRATION 
Multicomponent Mixture . The no-slip boundary condition for the species 
concentration without any assumptions may be obtained from Eqs. (44) and 
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(45). In the absence of slip, the Knudsen layer thickness shrinks to zero 
and the values at the top of the Knudsen layer become the wall values (see 
Fig. 1): 
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where we have neglected the higher order shear (i.e. P^/p * 1). 

It is suggested here, again, that the concentration for the major 
specie (for example, nitrogen) be obtained by requiring the sun of concen- 
trations of all the species to equal unity. 

For the recombining 0 and N atoms, Eq. (66) may further be simpli- 
fied to: 
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( 68 ) 


The recombination rate constant k ^ in Eq. (68) has been defined as (Ref. 

2 ) 
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f| = incident flux 

F* = specularly reflected flux 

F w = diffusely reflected flux 

r i 



Figure 1. The Knudsen layer showing general fluxes and coordinate axes. 

The temperature as a function of normal distance is schematically 
overlayed. 
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(69) 
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without slip and higher order shear. 

For a noncatalytic wall (with = 0) , Eq. (49) gives 

3 C 

( — = 0 (70) 

3y 

for all the species of a multicomponent mixture. 

Binary Mixture . For a two-species mixture of atoms and molecules, Eq . 
(55b) gives 


w 1 9 ^a 
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for a surface with finite catalycity and 
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for a noncatalytic surface. Eqs. (71) and (72) are employed in Ref. 4. 

Appendices D and E give the slip and no-slip boundary conditions (pre- 
sented in the text earlier), in the dimensionless form for the body-fitted 
and spherical polar coordinates, respectively, for a planar flow. The 
various integrals employed for evaluation the net, incident, and specularly 
reflected fluxes defined through Eqs. (5) to (7), respectively, are given in 
Appendix F. 
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DISCUSSION AND CONCLUSIONS 

The present analysis provides a consistent formulation for the slip 
equations for a multicomponent, binary, and single species mixture reacting 
catalytical ly on the surface. The differences between the slip expressions 
obtained by various researchers have been reconciled and implication of 
various assumptions (some of them inconsistent) contained in those express- 
ions is discussed. The slip equations have also been obtained in body-fit- 
ted and spherical polar coordinates in a form which can readily be employed 
in the flow-field calculations. 

Usually, the equation for pressure slip is not required as a boundary 
condition, but is needed to obtain the surface pressure. The temperature 
slip equation given here is for a constant surface temperature, which is 
provided as a boundary condition. For an adiabatic surface, however, the 
slip temperature, T , may be obtained by equating the wall heat-transfer 
rate to zero i.e. 


( „ 3T 

% = u — - 

3 y 


NS 
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i=l 


9v 


j i h i + W V 0x 


Ox 


3 y 


•) s = ° 


where the expression for mass flux j. is provided in Appendix B and the 
higher order terms have been dropped. The temperature slip equation will 
now be required to obtain the wall temperature, T . 

Tt 

An expression has also been obtained for the finite-rate species-con- 
centration boundary condition for a multicomponent gas mixture without sur- 
face slip. This boundary condition in the literature (Ref. 4) has generally 
been specified by assuming the dissociated air as a binary mixture. 
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However, the binary mixture formulation does not provide boundary conditions 
for the recombined molecules (on surface) such as O 2 and N 2 in a multicompo- 
nent mixture; it gives boundary conditions only for the recombining atoms. 

To assess the importance of various terms in the general boundary con- 
ditions suggested here, a study is being undertaken in which the flow prop- 
erties and boundary conditions would be varied systematically. The boundary 
equations form a simultaneous set, which is being coupled with a flow field 
calculation procedure in the stagnation region. These equations would 
finally be coupled with the viscous shock-layer code developed by Moss (Ref. 
13) for the detailed flow-field calculations. The boundary equations 
obtained in the present work should provide a more realistic set of boundary 
conditions for a multi -component mixture for low-Reynolds-number slip flows 
as well as no-slip flows. 

In conclusion, the boundary slip expressions obtained here are closed 
form solutions of the mass, momentum, and energy flux equations using the 
Chapman-Enskog velocity distribution function. This function represents a 
solution of the Boltzmann equation in the Navier-Stokes approximation. 
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APPENDIX A 


CONSTANTS APPEARING IN THE DISTRIBUTION FUNCTION 
AND DIFFUSION VELOCITY IN TERMS OF TRANSPORT PROPERTIES 

The Sonine expansion coefficients a. Q , a^, b. Q , and c!^ appearing 

in the general velocity distribution function are found by a variational 

technique in which they are solutions to sets of simultaneous equations. 

References 2 and 7 provide the coefficients in terms of solutions to these 

set of equations. These solutions are expressed in terms of the collision 

( s ) 

integrals &-jj* • Reference 7 also provides the transport properties in 

terms of the Sonine expansion coefficients. Thus, in place of evaluating 
these coefficients in terms of the collision integrals, they may be 
expressed in terms of the transport properties. The various relations 
are: 


( i ) for a 
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T n. m. f 2kT 

o{ (?) * - 



a i0< 5 > 
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( i i ) for a . 
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(iv) for c. 


( j) 
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In the distribution function used here the Kernel C^- C-^ 
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has associated Sonine polynomial coefficients c^’^. Now o (see 

Ref. 7), then c^ = c io’^ and 


'ij(^) 



(j) 


c i0 


(O 


(A4) 


where d| (? ) , K (K), m{K), and D-. (?) are the mul ti -component thermal 

diffusion coefficient, thermal conductivity, coefficient of viscosity, and 

the multicomponent diffusion coefficient, respectively. The argument ? is 

the number of terms used in the Sonine expansion. Except for dT, letting 

? = 1 gives quite good results for K, y and D... When ? = 1, however, 

^ J 

coefficients d| vanish. Hence, in order to get the coefficient of thermal 
diffusion, it is necessary to take at least two terms in the Sonine expan- 
sion (i.e. ? = 2). If the argument does not appear with a coefficient (ex- 
cept for a^), it is considered to have one term in the Sonine expansion. 

The diffusion velocity for ith species in terms of the transport co- 
efficients is obtained as (see ref. 7) 
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where the diffusion vector d^ for the jth species is defined after ne- 
glecting the external forces as (see Eq. (16) in the text): 


. ,, n . n . n .m . 

d J = (_i) + (^_ - — — ) (Jtnp) ( A6) 
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Eq. (A6) may be simplified to obtain 
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(A7) 


SIMPLIFIED EXPRESSIONS FOR COEFFICIENTS a^ AND b. Q 
Since a^ and b^ can not be obtained directly from equations (A2) 
and (A3), we introduce the following assumptions for the mixture thermal 
conductivity and viscosity; respectively. 


NS n. 

K - I 4 K. ( A8) 

i=l n 1 


NS n. 

v - i 4-^i ( A9 ) 

i=l n 


Equations (A8) and (A9) are approximate forms for the more exact formula of 
Wilke (see ref. 14). These equations imply that 

NS n . 

I 4 *ii - 1 (A10) 

j=l n J 

in the Wilke's relation. This is approximately true for air. 

Thus, with the help of Eqs. (A2), (A3), (A8), and (A9) we may obtain 
the following approximate expressions for a^ and b^, respectively. 


n 


i 




n. 

l 


n 



(All) 


43 



For accurate evaluations of a^ and b^, cumbersome expressions of the 
type of Eq. (7.4-49) of ref. 7 are required to be solved. 


APPENDIX B 


EVALUATED EXPRESSIONS FOR THE FLUXES OF MASS, 

MOMENTUM, AND ENERGY FROM INTEGRALS OVER 
THE DISTRIBUTION FUNCTION 

The fluxes of mass, momentum, and energy (i.e. of a property 
are given in terms of integrals over the velocity distribution function 
f(V) in the "Analysis" section of main text. Here we give evaluated forms 
for these fluxes obtained from integration over the distribution functions. 
The various integrals needed in these evaluations are provided in Appendix 
F. 

(i) Expressions for Mass Flux of Species 


Net 
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- ~~~ M t- 1^- (..T)J a, 0 ♦ n c(0> dj], ; 


2 m, 


3x, 


k = x, y, z 


(81) 


Incident 


4 


"1 "i 1 + ^i0 ( ^0x_ t ^0z_ 2 


z^r 


m 1 


6 3 X 3 Z 3 y 


+ i m. 
2 


(B2) 


Specularly Reflected 



(33) 


45 


Diffusely Reflected 
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( B4) 


It may be mentioned that M^. as obtained here is similar to the mass flux 
(jj) expression of Ref. 14. 

(ii) Correlation with Blottner's Expression (Ref. 15) for Mass Flux 

If the thermal diffusion term is neglected in Eq. (Bl) and 
substitutions are made for c^ and d J k from Eqs. (A4) and (A7), Eq . (Bl) 
may be simplified to obtain^" 
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where VF is the molecular (or atomic) weight of the ith species and is 
related to the mass of the ith species, m-, through the relation 
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Here k is the Boltzmann constant and R is the universal gas constant. 
Equation (85) may also be written as 
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or, with the introduction of the multicomponent Lewis nunber, L.., defined 

i J 


*It may be noted here that D-j -j = 0. See reference 7 for details. 
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as 


L 1J* P c pV ( ‘’ Pr V'‘> 


Eq. (B7a) may be rewritten as 
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where 
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In general , 
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L. . = L. . ( C., M., M., L..) 
1J lj l 1 J 1 3 


See Eqs. (7), (8a), and (8b) of reference 15. 
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If the binary Lewis nunbers are assuned to be the same for ail the 

species, then 


L ij = hd= L ij< c i- M i-V 


and Eq. (B8) for Ab^- becomes 


Ab ij = Le. 
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[ - L « + L '«. c -. ] 
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and the mass flux Eq. (B7e) now becomes 
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Further, if the binary Lewis numbers, are constant for ■>; ; chs 


species, the term 
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is zero in Eq. (BIO). Therefore, 
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Thus, Eq. (BIO) gives the mass flux due to concentration gradients only with 

Ab*^ and Le^ defined through Eqs. (B9) and (Bllb), respectively . 

★ 

To utilize Eqs. (B9) and (Bllb), L. . are still to be evaluated 

' J 

through complex matrix inversions (ref. 15). If we now make a little 

* . 

stronger assumption such that L i j = L 12 = Le (same for all the species), 

c * 

we obtain a much simpler expression for Ab^ and Eq . (Bllb) becomes 
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Similarly the expression for Ab^. 


may be obtained from Eq . (B9) as 


ab* = Le, - Le [!l + (1 - _1) (1-C,)] (813) 

•J 1 = u 
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3C i aC i 

It may be noted that for the case when = 0, the term Le. in Eq . 

3y 3y 

(BIO) vanishes and Eq. (B12b) is not required. Further, employing a 
constant value for the Lewis number, Le (= p Pr £ 12 /u), does not imply Le^ 
to have a constant value. 

It should be pointed out here that for the case when L.. is the same 
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for all the species we have used the condition 


NS 




i K ^ 1 * 0 


3 C 


J „ 
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= ★ 

(and not just Ab.^ = 0 as employed in ref. 15) to obtain Eq. (811b) for 
Le.. . The present formulation appears to provide a more consistent 
expression for Le^ . 

Finally, Eqs. (BIO), (B12b), and (B13) for a binary (two species) 
mixture give 
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j. = — — Le ( — 1) ; (Fick's law for mass flux) (B14) 
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(iii) Expressions For Normal Momentum Flux 
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are terms in the stress tensor for the ith species. 
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Incident and Specularly Reflected 
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(iv) Expressions For Tangential Moment un Flux 
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(v) Expressions For the Energy Flux 
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APPENDIX C 


SLIP BOUNDARY CONDITIONS GIVEN BY DAVIS FOR A BINARY MIXTURE 
In this appendix, the dimensional form of the slip boundary conditions 
obtained by Davis (Refs. 8 and 9) is provided for the purpose of comparison 
with the slip expressions obtained in the main text for a binary mixture. 

The dimensionless form of the surface slip conditions provided in refs. 
8 and 9 are 

Concentration slip (ref. 8)* 



Pressure slip (ref. 


9) 


tt 



(Cl) 


(C2) 


(C3) 


+ Here the dimensionless quantities are denoted with a 'bar'. See Appendix 
D for details of nondimensional i zation and the coordinate system. 

^For_the reasons explained in the main text, we have employed 'iT' in place 
of 'y 1 in this expression. 
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Temperature slip (ref. 8) 



With the help of the nondimensionalizing quantities given 
the following dimensional forms of Eqs. (Cl) through (C4) 
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in refs. 8 and 9, 
are obtained: 
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★ 

where n is the coordinate normal to surface (see Fig. D1 in Appendix D) 
and n and m are the nunber density and mass of a mixture molecule, 
respectively. The constants ai , bi , and ci are defined as: 


31 = 0 7 ( — } 

16 \ 2 9 
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16 \J2 0 ' 
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APPENDIX D 


SLIP CONDITIONS IN THE BODY-ORIENTED COORDINATE SYSTEM 

In this Appendix we obtain the slip conditions for an ax i symmetric body 

* * 

in the body-oriented coordinate (s , n ) system. The coordinate configura- 
tion for such a system is given in Fig. Dl. The equations of this Appendix 
contain simplifications similar to those given for a multicomponent mixture, 
following Eq. (38c) of the main text. 



Figure Dl. Coordinate Configuration For A Body (Syimetric About Its Axis) 

★ 

By employing the metric coefficients (see ref. 16) of h x = 1 + n < and 
h 2 = 1, equations (39) through (45) give the following slip equations in the 
body-oriented coordinate system. 
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Concentration slip 
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where (from Eqs. (35a) through (35c)) we define M. as: 
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For all other atoms and molecules 
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Temperature slip 
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If the internal energy (comprising of rotational and vibrational energies in 
this case) is considered frozen during reflection from the surface, the 
following form of the temperature slip equation is obtained (see Eq. (28) of 
Ref. 5): 
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where 
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In Eqs. (Dl), (D5), (06) and (D8a) the ratio — is defined as 
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and bg is related to viscosity, u, through the relation 


y = I nkT b f 


(DIO) 


which has been obtained from (A3) by assuming b^ to be the same (bg) for 
all the species. 

Equations (Dl) through (D8) may now be nondimensional ized by using the 
following relations (refs. 8 and 13): 
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Introducing the nondimensional quantities as defined here, the following 
equations are obtained from Eqs. (D-l) through (D-8): 

Concentration slip 
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(D12b) 
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Velocity Slip 
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(D 1 3) 
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(D15) 
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p diatomic w i 
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k ur A cf 


Or, with frozen internal energy during reflection from surface. 
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Pr = 


v 

K 


(Prandtl nunber) 


pC D D 12 

Le = — - (Lewis nunber) 

K 

U. 

M = — (Mach nunber) 

oo " 


and 




W T T , 
00 s ref 


00 \ w a T « 



(Recombination rate 
constant) 


For a noncatalytic surface (y a = = 0), eq. (D18a) through (D18c) give 

^ = 0 for all species. For this case, Eq. (017) becomes 
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(D20) 


Equation (D20), similar to Eq. (47b), gives 



= 0 


( 021 ) 


which may be employed as the boundary condition for a noncatalytic surface 
with a multicomponent gas mixture (with the binary assumption for diffusion 
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coefficients). For a fully catalytic surface, employing = 1 and using 
Eqs. (018a) through (018c) in Eq. (D17) would yield the appropriate 
concentration si ip. 

If the multicomponent diffusion coefficients are retained without the 
binary assumption the underlined terms in Eq. (D13) are to be replaced by 
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and the underlined term in Eq. (014) are to be replaced by 
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where the multicomponent Lewis nunber, L. ^ , is defined as 

p C„ D- . 


L ij 


P LI 


Simplifications For A Binary Mixture 

When all the species in a gas mixture can be considered as atoms and mole- 
cules only (see ref. 8), Eqs. (Dll) through (D18) maybe further simplified. 
With the assumption of Py/p “ 1, Eqs. (Oil) through (D18) for a binary mix- 
ture yield: 
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Concentration Slip 



(D22) 


Pressure Slip 
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(D23) 


where we have used the equation of state 


1 r T s T refj j- 
M 2 . T m T, « s 


( D24) 


and neglected the higher order shear terms as well as diffusion terms. 


Velocity slip 




K U 




1+n < 


(D25) 


where we have again employed the equation of state given earlier and neglec- 
ted higher order shear, conduction, and diffusion terms. In obtaining 
Eq. (D25) we have also used the approximation: 


69 
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1/2 

x , t C A S t W M ' V - V 
[ c a - aTO 


To be more exact one may keep this factor in Eq. (D25). 


Temperature slip 




p w ^wA 

Xi 





( cl + 1 ) 
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( D26) 


where 

rc A (w M - w A > ♦ w / 2 

xi * — 

3/2 3/2 3/2 

[C A (W M - 2 W A ) + 2 W fl ] 

and — < Xi < 1 for 0 < < 1 and diatomic molecules (W^ = 2 W^) . In 

obtaining Eq . (026), we have employed the equation of state and = 2W ft . 

With frozen internal energy during reflection from the surface, the 
temperature slip equation becomes 
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where 


[C fl (w„ - W A ) * 

r 3/2 3/2 3/2 

[ C A (W M - W A ) + W A ] 


For diatomic molecules (W» = 2 W A ) , X 2 may be taken as unity for 
0< C A < 1. Equation (D27) employs = 2W A> 


Equation for C A 



Equation (28), for a noncatalytic surface with y a = k wA = 0, gives 
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(D27) 
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(D29) 
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whereas using Y ^ = 1 in Eq. (D29) and in the expression for k : 



(D30) 


would yield the appropriate concentration slip for a fully catalytic sur- 
face. 

With the following (somewhat inconsistent) assumptions, Eqs. (D23), 
(D25) , (027), and (D28) maybe simplified to those obtained in refs. 8 and 
11 : 

(i) In pressure slip Eq. (023), Pr = 1 i s employed along with approxi- 
mations 


and 


i (JL), . 3 

5 7-1 


(V 2 [ c a"a + (V 2 » -M - .r s 

W A “m 


(ii) In velocity slip Eq. (025) the following assumption is made for the 
mixture molecular weights 


W 
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(iii) In temperature slip Eq. (D27) an assumption is made for the surface 
to be noncatalytic so that k w ^ = 0. In addition, Prandtl number 
(Pr) is taken unity and it is assumed that 


I 



15 

8 


(iv) No assumption is required in the concentration slip Eq. (D28). But 
it may not be appropriate to employ (Ref. 11) this equation for all 
the species of a multicomponent mixture. This equation is applic- 
able to concentration slip of recombining atoms only in a binary 
(two-species) mixture (Ref. 8) of atoms and molecules. 

No equation has been obtained in reference 8 or 11 to correspond to Eq . 

w 

(D22) to obtain wall values of the species concentration (C^) from th§ 
values at the edge of the Knudsen layer (C^) . 


Slip Expressions For A Single Species Mixture 
For a single species mixture (y^ = 0), the following slip expressions are 
obtained : 


Density Slip 



(D31) 
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Pressure SI ip 



( 032 ) 


Velocity Slip 


u 5 



K U 


1+n < 


(D33) 


Temperature Slip [For a gas consisting of molecules only (i.e. diatomic 
perfect gas)] 



where we have used xi = — » or with frozen internal energy during 

2 

reflection from the surface. 



(D35) 


Equations (D32) (D33) and (035) are the ones employed in References 8 and 11 
with Pr = 1. 

No-Slip Species Concentration Boundary Condition 
Multicomponent Mixture 

The no-slip boundary condition may be obtained from Eqs. (D17) and 
(D18) . in the absence of slip, the Knudsen layer thickness shrinks to 
almost zero, the values at the top of the Knudsen layer become the wall 
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values (See Fig. 1) : 
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_ xw . r A = 0 for M = O 2 
A ' l A » N for M = N 2 


^ = 0 for all other species 


(D37b) 

(D37c) 


where we have neglected the higher order shear (i.e. Py/p W «* 1). 

For the recombining 0 and N atoms, Eq. (D17) may also be written as 
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The recombination rate constant k wA in Eq . (D37a) and (D38) has been 
defined as (Ref. 2) 
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without slip and higher order shear. 

For a noncatalytic wall (with y a = k wA =0), Eq. (D21) gives 


3C 

(— L = o 


(D40) 


for all the species of a multicomponent mixture with the binary assumption 
about the diffusion coefficients. 

Binary Mixture 

For a two-species mixture of atoms and molecues, Eq . (D28) gives 



(D41a) 


1 


which may also be obtained from Eq. (038) for a surface with finite 
catalycity. k wA in Eq . (41a) is again obtained from Eq . (D39). For a 
noncatalytic surface with F wA = 0, Eq. (D41) gives 

3 C 

hr)„ = 0 <°«> 

3 n 
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Eqs. (D41) and (042) are obtained in ref. 4. However, ref. 4 employed Eq . 
(030) in place of Eq. (039) for obtaining the reaction rate constant, k^. 
As obtained in ref. 2, Eq. (D30) is more appropriate when slip and higher 
order shear are included in the flowfield analysis (e.g. under rarefied or 
low density conditions). 

If k ^ is substituted from Eq. (039) in Eq. (D41a), we may also 
obtain 
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T, Pr 




(— ) w (D41b) 

3rT 


where we have used the equation of state 

P W 1 , T w . 

— (D43) 

m2 _ = 

p “> y T W 

w 'oo * vy 


Eq. (D41b) compares with the corresponding equation of ref. 8, if one 
keeps in mind that the diffusion coefficient Fp (of Ref. 8) in the absence 
of slip is related to the present variables by the relation 


p w - ’ r Le 'i 

f d y w y* 


It may be seen from Eqs. (41a) or (41b) that the gradient (3C./3n) is 

n W 

governed by the ratio or /e 2 . Therefore, for surfaces which are 

almost noncatalytic, this ratio would be of the order of one for large 
values of the Reynolds number parameter (1/e 2 ). This would imply that a 
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surface may behave effectively as more catalytic under high density than 
under low density conditions. 


appendix e 
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SLIP CONDITIONS IN THE SPHERICAL COORDINATES 
For analyzing the flow in the stagnation region of a body, it is conven- 
ient to use the spherical polar coordinates (r, 4 >) for the two-dimensional 

flow. The following relations exist between the spherical polar and the 
• * * 

body-oriented (s , n ) coordinates over the spherical portion of a body 



Figure El. Coordinate Configuration 


* * 

r sin <(> = R + n cos 3 

r = r N + n* = r N ( 1 +n*< ) ; * = * 

' . n r N 

r N i> = s* ; B = 90 - * 
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or, in the non-dimensional form. 


r sin <j> = R + n cos 4> 1 


r = r M + n = l + n ; < = 1 
N 

<j)-s ; 3r = 3 n ; 3<j> = 3 s 


(E2) 


where we have used nose radius, r^, to nondimensional i ze all the distances 
(see Appendix D). 

Using the relations given by (E2) in Eqs. (Dll) through (D19), we can 
obtain the slip equations in the spherical polar coordinates for a multicom- 
ponent mixture. Once again, these equations employ the simplifications 
given for a multicomponent mixture, following Eq. (38c) of the main text. 
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Pressure Slip 
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Velocity Slip 
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Temperature -Slip 
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or, with frozen internal energy during reflection from surface. 
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■pS t-5 . r A = 0 for M = 0 2 

’•'M " *A ’ 1 A = N for M = N 2 


(ElOb) 
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= 0 for all other species 
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In obtaining these various equations, the following form of the equation of 
state has been employed: 


p s _ 1 ( T s T refj 

7 Y. L Is 


(E12) 


Similar to Appendix 0, Eqs . (E3) through (Ell) may further be simpli- 
fied for a binary mixture and for a single species mixture. 
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APPENDIX F 
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1 
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INTEGRALS REQUIRED FOR OBTAINING VARIOUS FLUXES 
In order to obtain the normal fluxes of mass, momentum, and energy, it 
is necessary to evaluate various integrals over the velocity space of the 
distribution function. These integrals involve terms that are various 
velocity moments of the distribution function. The integrals are provided 
for the net, incident, and specularly reflected fluxes. These integrals are 
basically the same as those provided in Reference 2 except for some 
corrections and additions. 


FI. Integrals Over the Entire Velocity Space (For Net Quantities) 
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F2. Integrals Over the Lower Half Velocity Space (For Incident Quantities) 
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F3. Integrals Over the Upper Half Velocity Space (For Specularly 
Reflected Quantities) 
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( -W y) appearing in the integrals emphasizes that the sign of the thermal 
velocity component normal to the surface, V (and consequently that of W ) 

, y y 

changes in the distribution function f upon reflection from the surface. 
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